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Abstract 

The eigenfunctions and eigenvalues of the linearized Boltzmann equation for inelastic hard 
spheres {d = 3) or disks (d = 2) corresponding to d + 2 hydrodynamic modes, are calculated 
in the long wavelength limit for a granular gas. The transport coefficients are identified and found 
to agree with those from the Chapman-Enskog solution. The dominance of hydrodynamic modes 
at long times and long wavelengths is studied via an exactly solvable kinetic model. A collisional 
continuum is bounded away from the hydrodynamic spectrum, assuring a hydrodynamic descrip- 
tion at long times. The bound is closely related to the power law decay of the velocity distribution 
in the reference homogeneous cooling state. 

PACS numbers: PACS numbers: 45.70.-n, 05.20.Dd, 51.10.+y 
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The simplest model for a granular gas at low density is given by the Boltzmann equa- 



tion for smooth, inelastic hard spheres or disks ^|. It provides an appropriate context in 
which to address a number of fundamental issues. Primary among these is the existence of 
a macroscopic fluid dynamics analogous to the Navier-Stokes description for real gases. The 
derivation of hydrodynamic equations from the inelastic Boltzmann equation and identifi- 
cation of expressions for the transport coefficients has been a problem of interest for two 
decades BUflQ- 

Accurate prediction of the transport coefficients as a function of the restitution coefficient 
0], and confirmation via Monte Carlo simulation Q], has been accomplished only within the 
last few years. However, the method used in these derivations (Chapman-Enskog) is formal 
and has been susceptible to debate. It is of interest therefore to look for a hydrodynamic de- 
scription in the simplest context of linear response to a small perturbation of the "universal" 
homogeneous cooling state (HCS). This response is provided by the linearized Boltzmann 
equation, and the existence of hydrodynamics depends on the spectrum of the operator. An- 
other issue is the isolation of the hydrodynamic spectrum from the remaining "microscopic" 
excitations, such that the hydrodynamic modes dominate on the long space and time scales. 
For elastic collisions, this has been proven for "hard" potentials (including hard spheres) . 

Here a calculation of the hydrodynamic spectrum from the linearized inelastic Boltzmann 
equation is reported. The transport coefficients are identified and shown to agree with 
those from the Chapman-Enskog method. Moreover, some support for the isolation of the 
hvdrodynamic spectrum . provded by a„ exact analysis of a mode. Bo.t.mann equation B 
12(1 . The model exhibits the dominant features of granular gases, including a non-trivial HCS 
with algebraic decay at large velocities, discussed extensively for related Maxwell models 
[lol |. It is shown that the pure point spectrum for elastic collisions develops a continuum 
for any degree of inelasticity, extending toward the hydrodynamic spectrum. However, it 
remains bounded away from the latter by an amount that is controlled by the power law 
of the decay of the velocity distribution. This relationship of the spectrum to the reference 
HCS is unexpected, but it is argued that it is more general than the limitations of the model. 

The inelastic nonlinear Boltzmann equation for the density / (r, v,t) of particles of mass 
m at position r with velocity v at time t has the form 

(9, + vV)/ = J[/,/], (1) 
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where J[f, f] is a bilinear functional of / Q]. There is no stationary solution to this equation 
for an isolated system. However, a homogeneous coohng solution (HCS) is assumed to exist 
with a scaling property similar to the Maxwellian, 

fncsi^, t) = nHVo'^{t)xo{c), c = . (2) 



Here uh is the homogeneous density and fo(^) = \/2T(t)/m is the thermal velocity defined 
in terms of the temperature T{t). The latter obeys the equation T~^dtT = —(o, where the 
cooling rate Co = ClXcXo] also is a specified bilinear functional. Both xo and Co must be 
determined self-consistently from the Boltzmann equation 

1 d 

2^°9c ' ^^^^^^^ ^ J[fHcs-, fucs]- (3) 



The existence of fn is supported by approximate polynomial expansions 



Carlo simulation 
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and by Monte 



To study the relaxation of small spatial perturbations of the HCS, A is defined by 
/ = fucs [1 + ^]) and terms up through linear order in A are retained in the Boltzmann 
equation. The resulting equation is then written in dimensionless form with the velocities 
scaled relative to VQ{t) and the time expressed in terms of s = dt'vo {t'), where z/o(t) is 



the average collision frequency, proportional to \jT{t). It is sufficient to consider a single 
Fourier component A(k) that obeys 

((9, - ik ■ c + £) A = 0, (4) 

>CA = -Xo ^ {^[Xo, XoA] + J [xoA, xo 

-ll-(cXoA)}. (5) 

The solution is sought in a Hilbert space defined by the scalar product 

(a, 6) = j dcxo{c)a*{c)b{c). (6) 

The relevant eigenvalue problem is 

{-ik ■ c + C) (f), = uji (k) (j),. (7) 

It is understood that the index i may be discrete or continuous. The hydrodynamic spec- 
trum can be defined precisely as follows. The linearized Boltzmann equation provides exact 



balance equations for the moments of A corresponding to the density, flow velocity, and 
temperature. The latter are the d + 2 hydrodjTiamic fields. The spectrum of these bal- 
ance equations can be calculated without approximation in the limit k = 0. For elastic 
collisions, there is a (rf + 2)-fold degenerate point at zero eigenvalue, corresponding to the 
conservation laws for binary collisions. Their perturbation at finite k defines the hydrody- 
namic modes more generally. For inelastic collisions, the spectrum is again d + 2 points, 
now at 0, Co/2, —Co/2, with the latter being d-fold degenerate. Again, the hydrodynamic 
modes more generally are defined as those solutions to ((Tj) that are continuously connected 
as functions of k to these^pecial A; = solutions. For elastic collisions, it is possible to prove 
analyticity about k = Q [8J, and it is assumed this is the case here as well. 

The first step is to show that C includes the spectrum of the balance equations at A; = 0, 
and to determine the corresponding eigenfunctions. For elastic collisions, the stationary state 
is given by J[/a/5 /a/] = 0, whose solution is the Maxwellian. Differentiating this stationary 
equation with respect to the hydrodynamic fields leads to the linearized Boltzmann operator 
acting on combinations of 1, c, (? independently being zero. These are the d+2 eigenfunctions 
with zero eigenvalue. A similar analysis works for the inelastic case, using the stationary 
condition (jS)). Straightforward calculations give 

^# = ^f#, ^ = l,..,rf + 2, (8) 

0(o)^{d + l + c.^,-rf-c.^,-^j, (9) 

-.{0,C/2,-C/2}. (10) 

The evaluation of the hydrodynamic spectrum at finite k is now a well-defined technical 
problem. In general, the eigenvalues depend on both k and the restitution coefficient a. 
It is known from the Navier-Stokes equations that the behavior for small k {a) is not 
uniform with respect to a {k), so that a perturbation expansion in one or the other does 
not provide the entire mode structure. Here, for illustration, only the results for small 
k are reported. The expansion is defined by cjj (fc) = uf'^ + ikuf'^ + k'^uf'^ + ■ ■ and 
(pi {k) = (pf^ + ik(j)i^ + k'^'pf^ + ■ ■ ■■ The leading terms are those of Eqs. © and (fTO|l . Since 
C is not self-adjoint, it is necessary to introduce a set of functions {'ipi} that are biorthogonal 
to {(pi). The leading terms of a similar expansion in k are found to be 

V^f -{l,^ + ^,k.c,e«.c|, (11) 



where {k, e^*)} are d pairwise orthogonal unit vectors. The results to second order in k are 
= 0, 

ujP = {4'\ k-c0f + 6,2(0 [xo, XoQi^f ^] , (12) 

0« = Q.(/:-..f)~'k-c0f), (13) 
# = -Q.(/:--n"k-c#. (14) 

The operators Qi = 1 — Vi are projections orthogonal to (pf^ with Via = (pf'^ «)• The 
constants cj) are transport coefficients which are identified from a similar calculation using 
the Navier-Stokes equations. For example, the shear viscosity rj is 

This agrees with the result obtained from the Chapman- Enskog method [13]. A more com- 
plete connection with the latter can be established as follows. Assuming that the hydro- 
dynamic spectrum dominates for long times and small k, the solution to the Boltzmann 
equation for small perturbations becomes 

d+2 



f fhcs (v) 



(16) 



1 + 5: 0, (k) 5y1 (k, s) 

Here Sy^ (k, s) are the hydrodynamic fields 

6y^ (k, s) = A (k, s)) = e'^^^'^ A (k, 0)) . (17) 

The connection between these fields and the macroscopic density, temperature, and fiow 
velocity is obtained for small k by using ipi ipf'^ in (fTTjl . The fields 6y^ (k, s) represent 
components of the initial distribution that obey the hydrodynamic equations for all times, 
whereas the density, temperature, and fiow velocity become equal to these only on the long 
space and time scales ^. Use of 0j {k) —>■ (pf^^ + ik(f)\^^ for small k in (fTBj) leads to a result 
;hat is in agreement with the Navier-Stokes approximation of the Chapman- Enskog method 

The existence of a hydrodynamic spectrum is relevant only if the corresponding modes 
dominate for long times and long wavelengths. A sufficient condition is that the magnitude of 
the hydrodynamic eigenvalues be smaller than that of all other parts of the spectrum. This is 
the case for elastic collisions, but the proof does not extend directly to — ik ■ c+C here and the 
question remains open in general In fact, the mode at Co/2 gives reason for concern since 



it might become comparable to the microscopic excitations at large inelasticity. To explore 
this point, it is instructive to consider a model for the Boltzmann equation that preserves the 
essential physics of a non-trivial HCS, exact balance equations, and hydrodynamic modes, 
as described above. The Boltzmann collision operator can be written 



J[f, f] = -Hf- g) 



(18) 



where vf is the "loss" contribution and vg is the "gain" contribution. Both the collision 
frequency v and g are specified positive functionals of /. 

The class of models to be considered here are characterized by the following two sim- 
plifications: 1) V is independent of v and, 2) v and g are functionals of / only through 
its moments with respect to l,v, . The first condition is similar to that of "Maxwell 
models" that have been introduced recently . The second condition means that v and g 
depend on the state of the system only through the hydrodynamic fields. The forms of these 
functionals are subjected to preserve the exact properties of the collision operator required 
by the balance equations, 

/ 1 \ / ^ \ 

^(r,v,t) = 



mv 



J 



n 

nu 



uTd (a^ + ^) ; 



(19) 



Here u is the flow velocity and \^ = 1 — C,/^. In principle, C, is the same functional of 
/ as obtained from the Boltzmann equation, although it can be eventually approximated 
similarly to v and g. I any case, Eq. (01) must be preserved. In the following, we will assume 
that both C and u scale as T^/^ in order to mimic hard sphere behavior. With the above 
class of models, it is obtained from Eq. Q 



dxP{x){Xx) '^g*{c/Xx) 



(20) 



where P{x) = px"*^^'^^-' with p = IvjC,. The dimensionless gain function g (^v,y^{r,t)^^ 
nu {Xvof g* [c/ X) has been introduced. Further investigation requires specification of 
g*{c/X). However, it is easily verified that moments of degree equal or greater than p 
are not finite if the corresponding moments of g* exist. This is consistent with a similar 
result for the Maxwell models where algebraic decay of Xo (c) is observed for large velocities, 
also due to a slower decay of the loss term than the gain term. 
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Linearization of the kinetic model around this HCS leads again to Eq. (0]), but with a 
simpler form for C It has a decomposition into the subspace of and its complement, 
L = VCV + QCQ, where V = J^T^i and Q = 1 — V. In the subspace of (pf^ the spectrum 
is just the d + 2 points of (fTUI) . 

d+2 

VCV = ujfVi. (21) 

i=l 

In the complementary subspace the operator is 

QCQ = vQ + Xo ■ (cxoQ) . (22) 

Interestingly, there is no explicit dependence on the gain function g* except through xo- For 
elastic collisions, the second term is missing and the model becomes the linearized Bhatnager- 
Gross-Krook model Q|, with a single degenerate point in the spectrum representing all the 
microscopic excitations of the Boltzmann operator. It is somewhat easier to analyze the 
corresponding adjoint problem 

"""^ " ■ = = ^^^^ 

where Q^ is the adjoint of Q. Solutions to this equation exist from the class of homogeneous 
functions 

0(xncf, E?i = ?- (24) 

i=l i=l 

The property = can be satisfied with a suitable choice of the linear combinations. The 
corresponding eigenvalue is u; = z/ — gCo/2. Thus the spectrum of C includes a continuum 
with 3?C(j < V. Of critical interest for hydrodynamics is whether the lower bound for this 
continuum on the real axis intersects the hydrodynamic point spectrum. The bound is 
determined by the maximum degree of such homogeneous eigenfunctions admitted in the 
Hilbert space. The condition is that exists. Using Eq. (j2Uj) . this becomes 



/oo 
^^^2Kg-(l+p)^ (25) 

/•oo 

Mq = / dcc^^'^g* (c) . (26) 



For finite Mg, the value of q is bounded from above by 3?g < Qm = p/2 = z^/Co- The 
continuous spectrum therefore is restricted by z//2 < ^uj < v. This means that the discrete 
spectrum is isolated from the continuum for all Co/^ < 1- But C,q/v must be smaller than 
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unity because of the third condition required by Eq. (fT^ . This condition is related with 
the fact that only a fraction of the energy of the particles can be lost in collisions. The 
continuum is singular in the sense that it is independent of the restitution coefficient but it 
is present only for finite inelasticity. 

For this class of models, it is seen that the hydrodynamic spectrum at /c = is isolated 
from the rest of the spectrum. The model also allows exact calculation of the hydrodynamic 
eigenvalues Ui {k) at finite k, without the restrictions of the perturbation theory calculation 
above. For elastic collisions, such a calculation shows that the hydrodynamic modes exist 
only for k less than some maximum value. It is expected that a similar condition applies in 
the granular well. 
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